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Abstract. Using a sample of ten nearby (z < 0.15), relaxed galaxy clusters in the temperature range [2 - 9] keV, we have 
investigated the scaling relation between the mass at various density contrasts (8 = 2500, 1000, 500, 200) and the cluster 
temperature. The masses are derived from NFW-type model fits to mass profiles, obtained under the hydrostatic assumption 
using precise measurements, with XMM-Newton, at least down to S = 1000. The logarithmic slope of the M-T relation is well 
constrained and is the same at all <5, reflecting the self-similarity of the mass profiles. At 6 = 500, the slope of the relation for 
the sub-sample of hot clusters (kT > 3.5 keV) is consistent with the standard self-similar expectation: a = 1.49 ± 0.15. The 
relation steepens when the whole sample is considered: a = 1.71 ± 0.09. The normalisation of the relation is discrepant (by 
~ 30 per cent), at all density contrasts, with the prediction from purely gravitation based models. Models that take into account 
radiative cooling and galaxy feedback are generally in better agreement with our data. We argue that remaining discrepancies, 
in particular at low 6, are more likely due to problems with models of the ICM thermal structure rather than to an incorrect 
estimate of the mass from X-ray data. 

Key words. Cosmology: observations, Cosmology: dark matter, Galaxies: cluster: general, (Galaxies) Intergalactic medium, 
X-rays: galaxies: clusters 
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From a theoretical point of view, galaxy clusters are charac- 
terised by their mass. Models of structure formation predict 
the space density, spatial distribution and physical properties of 
clusters (internal structure, radius, tem perature, luminosity , etc) 
as a function of mass and redshift (see Bertschinger 1998, for a 
review). However, observationally, the mass is not easily mea- 
sured, and the observed scaling relations are in fact expressed 
in terms of the temperature T, rather than the mass M. These 
scaling relations are important so urces of information on the 
physics of cluster formation (e.g. lVoit & Ponman1l2003l) . For 
the information to be complete, we must determine the M-T 
relation itself, which provides the missing link between the gas 
properties and the mass. Furthermore, measures of the cosmo- 
logical parameters, such as erg, Q m and w, from cluster abun- 
dance or spatial distribution, rely heavily on this relation to link 
the mass to the X-ray observables available from X-ray cluster 
surveys. The present error on the value of <x 8 , as determined 
from X-ray observatio ns, is dominated by uncertainty on the 
M-T relation JPierpaoli et alJl2003l Iviana et alJl2003l iHenrvl 
l2004h . and a precise calibration of this relation is ma ndatory if 
we w ant to do 'precision' cosmology with clusters JBor gani 
l2003h . 



Send offprint requests to: 
marnaudOdiscovery . saclay .cea.fr 



M. 



Arnaud, 



e-mail: 



The average temperature is expected to be closely related to 
the mass, via the virial theorem. We can first define M s as the 
mass within the radius R s , inside which the mean mass den- 
sity is 6 times the critical density, p c (z) = 3/z(z) 2 Ho/8^G, at 
the cluster redshift. We then expect h(z)Mg = A{6)T 3 ^ 2 , if clus- 
ters are in hydrostatic equilibrium and they obey self-similarity. 
Here h(z) is the Hubble constant normalised to its local value 
and A{5) depends on the internal structure. The above relation 
is remarkably well verified by adiabatic numerical simulations, 
down to 6 ~ 200, w hich roughly correspo nds to the virialised 
part of clusters (e.g. lEvrard & Gioial2002l) . 

For relaxed clusters, the mass can be derived from X-ray 
observations of the gas density and temperature profile and the 
hydrostatic equilibrium equation. In recent years, a sustained 
observational effort to measure the local M-T relation has been 
undertaken using ROSAT, ASCA and BeppoSAX, but no defini- 
tive picture has yet emerged. It is unclear whether the mass 
scales as T 3 ^ 2 as expected llHorne r9talJ fi999t lEttori et all 
120021 ICastillo-Morales & Schindlerl 1200311 : or if this is true 
only in the high mass regi me (kr > 4 keV), with a steep- 
ening at lower mass jNevalainen et"alJl200(j Ixu et alJl200ll 
iFinoguenov et alJl200ll): or even if the sl ope is higher over the 
entire mass range JSanderson et al.l2003l) . The derived normal- 
isations of the M-T relation derived from ASCA data are gen- 
erall y lower than predicted by adiabatic numerical s imulations 
(e.g. Nevalainen et al . 2000; Finoguenov et al .1200 lh . typically 
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40 per cent below the prediction oflEvrardetalJ i 19961) . On the 
other hand, using BeppoSAX data. lEttori*eTa?T l l2002l) found a 
normalisation consistent with the predictions (although the er- 
rors were large). 

These studies had to rely largely on extrapolation to derive 
the virial mass, and were limited by the low resolution and sta- 
tistical quality of the temperature profiles. With XMM-Newton 
and Chandra we can now measure the mass profile of clus- 
ters with unprece dented accuracy. Using Chandra observa- 
tions, lAllen et ail d200 lh derived an M-T relation slope of 
1.51 + 0.27, consistent with the self-similar model, and con- 
firmed the offset in normalisation. However, their sample com- 
prised only 5 hot (i.e massive) clusters (kT > 5.5 keV), and, 
due to the relatively small Chandra field of view, their M-T 
relation was established at Ry.^oo (i.e. about ~ 0.3j??.no). 

In a recent paper ( Pointec outeau. Arnaud & Pratn 12005 . 
hereafter, Paper I), we measured the integrated mass profiles of 
ten relaxed, nearby clusters observed with XMM-Newton. The 
sample has an excellent temperature coverage, from 2 to 9 keV. 
The mass profiles cover a wide range of radii (from O.OI.R200 to 
0.72?20o)> and are particularly well constrained between O.l/^oo 
and 0.5/?200- In Paper I, we studied the structural properties of 
the mass profiles, in order to test current scenarios for the Dark 
Matter clustering. In this paper, these data are used to estab- 
lish a precise M-T relation up to the virial radius. In Sect. [3] 
we describe how we derive the temperature and mass data. In 
Sect. |5] we present and compare the M-T relations at various 
density contrasts. We discuss the reliability of the X-ray mass 
estimates in Sect. 0]. The derived M-T relations are discussed 
with respect to pie-Chandra/XMM-Newton results in Sect. [5] 
and with expectations from models in Sect. [6] Our conclusions 
are presented in Sect. [7] 

Throughout the paper, results are given for the currently- 
favoured ACDM cosmology, with Hq = 70 km/s/Mpc, Q,„ = 
0.3 and Q A = 0.7. 

2. The data 

2.1. The mass at various density contrasts 

We used the mass profiles determined in Paper I to estimate the 
mass at four density contrasts: 6 = 2500, 1000,500 and 200 1 . 
We recall that the mass profiles were derived from the observed 
density and temperature profiles (corrected for PSF and pro- 
jection effects) using the hydrostatic equilibrium equation. The 
mass and errors at each radius of the temperature profile were 
calculated using a Monte Carlo method. 

In Paper I, we found that t he mass profiles are well de- 
scribed by an NFW-type model jNavarro et all 19971) : 

ln(l + c 200 x) - c 2 q Q x/(1 + c 2 oo*) 



M NFW (r) = M : 



200" 



(1) 



ln(l + c 2 oo) - c 2 oo/(l + c 2 oo) 
where x = r//?2oo- The concentration C200, and the mass M200, 
are free parameters determined by the fitting procedure. /?2oo 

1 We chose these <5s for the following reasons: 6 = 2500 allows 
a direct comparison with previous Chandra results; 6 = 1000 is the 
density contrast limit of our observations; 6 = 500 corresponds to 
the edge of the virialized part of clusters in a conservative approach; 
6 = 200 is the classical 'virial' radius in an £1 = 1 universe. 



is related to M 20 o via M 20 o/^2oo = (47r/3)200p c (z). The best 
fitting M200 and /?2oo values are recalled in Table ^ The con- 
centration parameters are given in Paper I. The correspond- 
ing NFW mass , Ms, at any density contrast 6, can be de- 
rived from the best fitting C200 and M200 values using Eq. ^ 
(M 6 = M N fw(^) with Mslr\ = (4n /3)6p c (z)). However, the 
computation of the error is not straightforward since the uncer- 
tainties in C200 and M200 are correlated. To avoid this problem, 
for each 5 under consideration, we refitted the observed mass 
profile data with an NFW model using Ms and the correspond- 
ing concentration as free parameters in the fit. The errors on Mg 
can then be derived from standard A^ 2 criteria for one inter- 
esting parameter, the concentration being optimised (minimum 
X 1 ) for any given value of Ms. The resulting M2500, A^iooo and 
M500 values are also shown in Tabled We also list the outer- 
most radius (in units of R200X and the corresponding density 
contrast (5 0DS , reached by the temperature profiles (and thus the 
outermost extent of the measured mass profiles). 

We will use these mass estimates to study the Ms-T rela- 
tion at various density contrasts. Our study is thus based on 
a parametric model of the observed mass profiles, rather than 
directly on the measured mass data. Let us discuss this point 
in more detail. All of the clusters are observed down to at 
least 5 = 1000 (#1000 = (0.47 + 0.02) /? 2 oo, averaged over the 
whole sample), the only exceptions being A 198 3 and MKW9 
(5 obs ~ 1400). At 8 = 1000 and 5 = 2500, using the best fit- 
ting model rather than the data is simply equivalent to 'smooth- 
ing' the data (without data extrapolation). We checked that this 
does not introduce a bias in the following way: for each clus- 
ter, we estimated the mass at 8 — 2500 by interpolating the 
observed profile expressed as a function of density contrast 2 in 
the log-log plane. We then compared the interpolated value to 
that derived from the NFW fit to the mass profile. In all cases 
the values are consistent within their lcr errors. The ratio of the 
two values has a median value of 0.99 across the sample, and 
there is no significant correlation with mass. This reflects the 
fact that the NFW model is a good fit to these data, particularly 
in the 0.1/?20o - 0.5/?200 range (see Paper I). 

As explained above, the estimates of M2500 and M1000 are 
made (almost) without data extrapolation. However, the mass 
estimates at 6 = 500 and 8 = 200 do involve extrapolation of 
the data. The M500 and M200 estimates rely on the assumption 
that the best fitting NFW model remains a good representation 
of the cluster mass profile beyond <S b s . We further discuss the 
reliability of this assumption in Sect.0] 

2.2. Overall temperature 

To investigate the M-T relation, we need to define a global 
temperature. For this quantity, we used the overall spectro- 
scopic temperature of the 0.1/? 200 < r < 0.5/? 200 region. The 
lower boundary of 0.1 R200 was chosen so as to avoid most of 
the cooling core, where a large dispersion is observed in the 
temperature profiles (Fig. The upper boundary is limited 
by the quality of the spectroscopic data. An upper boundary 



2 M(6) is readily derived from M(r), using 5 = 'iM{r)IAnp c (z)r i at 
each discrete radial value of the mass profile. 
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Table 1. Physical cluster parameters. Masses are in units of 1O 14 M , and are given for a ACDM cosmology with Q m = 0.3, 
£2a = 0.7, Hq = 70 km/s/Mpc. Errors are lcr errors. 



Cluster 


z 


kT (keV) 


R 2 ao (kpc) 


M 2 500 


Miooo 


M 500 


M 2 ao 


Robs/RlOO 


Sobs 


A1983 


0.0442 


2.18 ±0.09 


1103 ± 136 


0.43 


±0.09 


0.77 ± 


0.22 


1.09 


±0.37 


1.59 ±0.61 


0.38 


1455 


MKW9 


0.0382 


2.43 ± 0.24 


1006 ± 84 


0.41 


±0.07 


0.66 ± 


0.14 


0.88 


±0.20 


1.20 ±0.30 


0.41 


1401 


A2717 


0.0498 


2.56 ± 0.06 


1096 ± 44 


0.45 


±0.04 


0.79 ± 


0.08 


1.10 


±0.12 


1.57 ±0.19 


0.54 


727 


A1991 


0.0586 


2.71 ±0.07 


1106 ±41 


0.58 


±0.05 


0.91 ± 


0.09 


1.20 


±0.12 


1.63 ±0.18 


0.60 


655 


A2597 


0.0852 


3.67 ± 0.09 


1344 ± 49 


1.08 


±0.07 


1.69 ± 


0.14 


2.22 


±0.22 


3.00 ± 0.33 


0.57 


713 


A 1068 


0.1375 


4.67 ±0.11 


1635 ± 47 


1.47 


±0.07 


2.69 ± 


0.16 


3.87 


±0.28 


5.68 ± 0.49 


0.58 


622 


A1413 


0.1430 


6.62 ±0.14 


1707 ± 57 


2.33 


±0.13 


3.66 ± 


0.27 


4.82 


±0.42 


6.50 ± 0.65 


0.79 


339 


A478 


0.0881 


7.05 ±0.12 


2060 ± 112 


3.12 


±0.31 


5.43 ± 


0.70 


7.57 


± 1.11 


10.8 ± 1.8 


0.58 


650 


PKS 0745-191 


0.1028 


7.97 ± 0.28 


1999 ± 77 


3.32 


±0.23 


5.41 ± 


0.49 


7.27 


±0.75 


10.0 ± 1.2 


0.57 


694 


A2204 


0.1523 


8.26 ± 0.22 


2075 ± 77 


3.62 


±0.22 


6.11 ± 


0.51 


8.39 


±0.81 


11.8 ± 1.3 


0.61 


580 



Columns: (1) Cluster name; (2) Redshift; (3) Spectroscopic temperature of the 0.1 < r < 0.5 R 2 qq region; (4) R 2 oo, roughly the virial radius in 
numerical simulations; (5,6,7,8) Total mass at density contrast 6 = 2500, 1000, 500, 200, derived from an NFW fit to the observed mass profile; 
(9) Fraction of R 200 spectroscopically observed (outer radius of the final temperature bin); (10) Corresponding density contrast. 
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Fig. 1. Temperature profiles. The temperatures have been nor- 
malised to the spectroscopic temperature measured in 0.1 < 
r < 0.5 R200 region; the radius has been scaled to /?200- The 
profiles have been corrected for PSF and projection effects (see 
Paper I for details). 



of 0.5 /?200 appeared a good compromise. Only the data from 
A1983 and MKW9 do not quite reach this radius; they are how- 
ever detected up to ~ 0.4 R200 (see TableQ. Note that 0.5 R200 
corresponds roughly to 5 — 1000. 

For each cluster, we performed an isothermal fit of the spec- 
trum extracted within the [0.1 - 0.5]/? 200 range, /?2oo being de- 
rived from the best fitting NFW model (see above). In the fit, 
the abundance was let free and t he N» was fixed to t he 21 cm 
value (except for A478, see lPointecouteau et alJ2004l) . We cor- 
rected the derived value for PSF blurring and projection effects 
using the ratio of the mean emission-measure weighted value 
of the temperature profile in the [0.1 - 0.5] /?2oo region after 
PSF/projection correction to the mean value before correction 
(see Paper I for details on the correction procedure). The cor- 
rection factor is generally negligible and is always less than 5 
per cent. The resulting temperature values are given in Tabled 



We could have estimated a 'mass-weighted' temperature in 
the 0.1/?200 < r < 0.5/?2oo region from the temperature profile. 
However, this temperature would still be a 'spectroscopic' tem- 
perature since it would be derived from averaging over mea- 
sured X-ray temperatures. It would not, strictly speaking, be 
equivalent to the 'mass-weighted' temperature derived from 
numerical simulations. We thus preferred to use the overall 
spectroscopic temperature, which is a directly measured quan- 
tity and can also easily be estimated in numerical simulations. 
Since the region is defined in scaled radius, it can be derived 
from the simulated temperature profil es, using for ins tance the 
approach proposed recently by Mazzo tta et alJ (120041) . In addi- 
tion, we note that only a global spectroscopic temperature can 
be usually estimated for high z clusters. In our approach, the ex- 
traction region can be similarly defined and our definition thus 
allows a consistent study of the evolution of the M—T relation. 

We note that lAllen et alJ d200ll) use a mass-weighted tem- 
perature, r25oo> estimated from the temperature profile in the 
r < /?2500 ~ 0.3/?200 region. In practice, their definition is 
equivalent to ours because i) the temperature profiles are fairly 
flat beyond the cooling core region (r > 0.1/? 200) in both stud- 
ies, and ii) the cooling core does not contribute muc h in mass 
to the average. This can be checked from Fig. 1 of lAllen et alJ 
(2001), where ^300 cannot be distinguished from the spectro- 
scopic temperature of the region beyond the cooling core. 



3. The M-T relation 

3.1. The M2500 - T relation 

In order to check the consis tency of ou r XMM -Newton results 
with the Chandra study of lAllen et alJ d200 lh . we first inves- 
tigated the M2500 - T relation. For our sample, 6 = 2500 cor- 
responds to an average radius of [0.29 + 0.02] Rioo, where the 
mass is particularly well constrained for all clusters. 

The sample studied by lAllen et all fcOOll) comprises hot 
lensing clusters (5.5 to 15 keV). We thus considered only the 
sub-sample of clusters with moderate to high temperatures (i.e. 
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Table 2. Results of power law fits to the M s -T and Rg-T relation at various density contrasts 6. The data are fitted with a power 
law of the form h(z)M s — A s x (kT/5 keV) ff and h(z)Rs = B 5 x (kT/5 keV*f, where kT is the overall spectroscopic temperature 
of the [0.17? 200 _ 0.5i?2oo] region. A ACDM cosmology is assumed: Q m = 0.3, Q.\ = 0.7, and Hq = 70 km/s/Mpc. 







M s - 


-T relation 








Rg-T relation 


6 




a 




Cfogjnt 


^(dof) 


nhp 


B s (kpc) 


s 


Whole sample 


















200 


5.34+ 0.22 


1.72± 0.10 


0.077 


0.051 


14.49(8) 


0.07 


1674±23 


0.57+0.02 


500 


3.84± 0.14 


1.7 1± 0.09 


0.064 


0.039 


12.65(8) 


0.12 


1104+13 


0.57±0.02 


1000 


2.82± 0.09 


1.71± 0.08 


0.053 


0.027 


10.74(8) 


0.22 


791±8 


0.57+0.02 


2500 


1.69+ 0.05 


1.70± 0.07 


0.041 


0.016 


9.33(8) 


0.32 


491 ±4 


0.56+0.02 


T > 3.5 keV 


















200 


5.74+ 0.30 


1.49± 0.17 


0.081 


0.064 


10.45(4) 


0.03 


1714±30 


0.50+0.05 


500 


4.10± 0.19 


1.49± 0.15 


0.064 


0.046 


8.34(4) 


0.08 


1129+17 


0.50±0.05 


1000 


3.00± 0.12 


1.49± 0.14 


0.048 


0.027 


5.91(4) 


0.21 


807+10 


0.50+0.04 


2500 


1.79+ 0.06 


1.51± 0.11 


0.025 




2.50(4) 


0.65 


500±5 


0.50+0.03 



Columns: (1) Density contrast S; (2, 3) Intercept and slope for the Mg-T relation: h(z)Mg = Ag x (kr/5keV)" with standard errors; (4, 5) 
Raw and intrinsic scatter about the best fitting relations in the log-log plane (see Sec. l3.2t ; (6) Chi-squared and degree of freedom; (7) Null 
hypothesis probability associated with the best fit; (8,9) Intercept and slope for the Rg-T relation: h(z)Rg = Bg x (kr/SkeVV 3 ; 



T > 3.5 keV), and fitted the M25Q0-T relation using a power 
law model of the form: 



Kz)M s 



kT 



5keV 



(2) 



Here and in the following, the fit is performed using linear re- 
gression in the log-log plane, and the goodness of fit is calcu- 
lated using ax 2 estimator taking into account the errors on both 
mass and temperat ure. We us ed the routine FITEXY from nu- 
m erical recipe s JPress et al. Ill992h . Note that, as in the study 
of lAllen et al. ( 1200 II) . the masses are scaled by h(z), which 
corrects for the evolution expected in the standard self- similar 
model. This scaling factor is small in our z range (h(z) = 1.07 
at z - 0.15) but varies between ~ 1.05 and ~ 1.28 for the 
Chandra clusters located at higher redshifts (0.1 < z < 0.46). 

The data are well fitted by a power law (^- 2 /dof = 2.5/4). 
The slope, a — 1.51 + 0.11, is perfectly consistent with the 
expectation from the standard self-similar model, and with the 
results from Chandra observations (a = 1.51 + 0.27). The de- 
rived normalisation, A = (1.79 + 0.06) x 10 14 M , is also con- 
sistent with the Chand ra normalisation (see Table[3J. As noted 
hv lAllen et all (kOOlh . such a normalisation is discrepant with 
the value derived from numerical simulations including grav- 
itational heating only: our m easured value is ab out ~ 30 per 
cent below the prediction of Evrard et all Jl996l). When the 
Chand ra data for 4 of the 5 clusters studied by lAllen et al] 
(2001) are added to the present data set 3 , the best fitting val- 
ues are almost unchanged (a = 1.52 + 0.1 with the same inter- 
cept). This is due to the larger uncertainties in the Chandra tem- 
perature and mass determinations compared to those measured 
here (see Fig. [2). Figure [2] shows the best fit for the combined 
XMM-Newton and Chandra data compared to the expecta- 
tions f rom the adiabatic numerical simulations of lEvrard et alJ 
J1996I) . 



3 The fifth cluster is PKS 0745-191, which is common to both sam- 
ples. We use only the XMM-Newton measurement here. 



,15 



5=2500 




kT (keV) 

Fig. 2. The M- T relation at 6 = 2500 as seen by XMM-Newton 
from the observation of 6 hot (kT > 3.5 keV), relaxed clus- 
ters. Filled squares show the XMM-Newton data points; the full 
line shows the best fitting power law. The data on 4 published 
Chandra clusters (triangles) have been added to the fit but due 
to their larger uncertainties, they do not change the parameters 
of the fit to the XMM-Newton data only (see text). The dashed 
black line is the pr ediction from adiabatic numerical simula- 
tions JEvrardetal.il 9 96). 



Still working at 5 — 2500, we performed a fit over the 
whole XMM-Newton sample, i.e. now including the four low 
mass systems. We obtain a = 1.70 + 0.07, and a normalisation 
A = (1.79 ± 0.06) x 10 14 Mq. The fit is acceptable, although 
formally less good (x 2 /dof = 9.33/8). The slope now differs 
significantly from the expected value of a = 1.5, and is just 
barely consistent with it at a 3<r level. This is further discussed 
in Sec. 16. 21 
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Fig. 3. The M - T relation as seen by XMM-Newton from a sample of 10 clusters covering a temperature range from 2 to 9 keV. 
From top to bottom and left to right, the M - T relation is given at the density contrasts 6 of 200, 500, 1000 and 2500 with 
respect to the critical density of the Universe. Measurements are plotted with error bars. In each panel, the best fit for the whole 
sample is overplotted as a solid line, and the best fit for the hot cluster subsample is plotted as a dotted line. The predicted relation 
from adiabatic numerical simulations JEvrard et al Jl99(j) is overplotted as a dashed line. The long-dashed line (panels 5 = 500 
and 5 = 2500) is the relation derived from a numerical simulation including radiative cooling, star formation and SN feedback 
jBorgani et all2004 . 



3.2. The M s -T relations up to the virial radius 

Figure shows the Mg-T relations at various 6, together with 
the best fitting power law (Eq. |2j in each case, and t he pre - 
diction from the numerical simulations of lEvrard et alJ {l996). 
The best fitting slopes and normalisations are listed in Table|2] 
together with the standard errors. The best fits are listed and 
plotted both for the whole sample, and for the sub-sample of 
hot clusters. The corresponding Rg-T relations are also given 
in the Table. 

The normalisation and slope are nearly independent param- 
eters for the whole sample. The covariance in log{Ag) and a, 
normalised to the product of their standard errors, is small: 
0.045, 0.071, 0.092 and 0.1 1 for 5 = 2500, 1000, 500 and 200, 
respectively. This is due to our choice of the pivot of the Mg-T 
relation, kT = 5 keV (Eq.[2}, close to the mean temperature of 
the whole sample (4.8keV) or the median value of its tempera- 



ture range (5.2keV). The normalisation for other pivots can be 
derived using our best fitting values from Eq. |5] For instance, 
for a pivot at 6 keV, used in several works on cluster scaling 
relations, the normalisation Ag is ~ 36 per cent higher. The rel- 
ative error, cr^JAg, is increased by ~ 10 per cent. 

We have also computed the raw and intrinsic scatter about 
the best fitting relations in the log-log plane. They are given 
in Table |2] To estimate the raw scatter, we used the orthogo- 
nal distances to the regression line, weighted by the error 4 . 
The intrinsic scatter is computed from the quadratic difference 
between the raw scatter and the scatter expected from the sta- 
tistical errors. 



4 For a linear relation of the form Y = aX + b, and a sample of N 
data points (K,, X,) with errors cry. and a Xi , an estimate of the square 
of the raw scatter is: cr 2 aw = -J^ Y*1Li w i(Yi ~ a X - b) 2 , where w, = 

( , /w) g; 1/gf with a) = cr\ + a 2 a%. Here Y = \og(M s ),X = \og(T). 
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The behaviour of the M2500-T relation is reproduced at all 
other density contrasts. The slope is stable on all spatial scales: 
the variation is at most 10 per cent of the statistical error. It is 
always consistent with the expected a = 1.5 value for the sub- 
sample of T > 3.5 keV clusters, whereas it steepens to a — 1.7 
when the cool clusters are included. Similarly, the normalisa- 
tion remains ~ 30 per cent below the value from adiabatic nu- 
merical simulations of lEvrard et al.l dl996l) at all 6. 

This stable behaviour is a direct consequence of the self- 
similarity of the mass profiles (Paper I). For an NFW type 
profile, the ratio of the masses at different density contrasts 
only depends on the concentration parameter, and Mg = 
M25mF(c,6). If all clusters had exactly the same concentra- 
tion parameter (i.e., perfect self-similarity) the Mg-T rela- 
tions at various 5 should differ only in their normalisation: 
M s = A(2500)F(c, 6)T a , for M 25 oo = A(2500)r ff . The clus- 
ters in our sample are not perfectly self-similar but there is 
no significant variation of c with mass, and thus temperature 
(see Fig. 3 of Paper I). This explains the observed invari- 
ance of the slope. Furthermore, the observed concentration pa- 
rameter is consistent with theoretical expectations (Paper I). 
As a consequence, the variation with 6 of the normalisation, 
A(6) = A(25QQ) F{c,S), follows expectations, and the offset 
with respect to simulations observed at S = 2 500 remains the 
same at all 5 (see also Fig. 0] and Fig 13 of IPratt & Ar naud 
120021) . 

However, the quality of the power law fit decreases with 
decreasing 6 (see Table |3 and also Fig. [3}- The reduced^ 2 in- 
creases and the corresponding null hypothesis probability for 
the whole sample varies from 0.32 at 6 — 2500 to 0.07 at 
6 = 200. This behaviour corresponds to an increase of the 
intrinsic scatter in the observed Mg-T relation (see Table |2}. 
For the whole sample, the worst fit and largest intrinsic scat- 
ter is observed at 6 — 200. The regression method we used 
(Sect. |2~T1 is strictly valid only if the intrinsic scatter is negli- 
gible as compared to the statistical scatter. This is not always 
the case (see Table[2), and we first checked if our results could 
be affected, using the M200-T relation for the whole sample 
(i.e., the worst c ase). We refitted the data us ing the orthogonal 
BCES method (Akrit as & Bershadvl fl996). While this is the 
least-biased regression method when both measurement errors 
and intrinsic scatter are present, it is less accurate that the \ 2 
method when the intrinsic scatter is negligible. The best fitting 
values remain unchanged (within 0.5 per cent) and the standard 
error estimates are only slightly larger (by ~ 15 per cent). 

The regression method we used is thus justified. 
Nevertheless, the derived intrinsic scatter should not be over- 
interpreted. The cluster sample is small and is certainly not 
representative of the entire cluster population. In particular, it 
is heavily biased towards the more relaxed clusters. Moreover, 
the increased intrinsic scatter at low 6 may be an artifact of 
the method we used to derive the various Mg. These were de- 
rived from an NFW fit to the observed mass profile, using Mg 
as a free parameter (see Sect. 12. 1> . In the NFW fit, M2500 is ex- 
tremely well constrained by the data around 5 = 2500, quasi- 
independently of the shape of the observed mass profile. On 
the other hand, Mg at low 5, beyond the maximum radius of 
observation, can be viewed as an 'extrapolation' of the NFW 



model best fitting the observed mass profile. It thus depends 
both on the normalisation of the observed profile (basically 
A/2500) an d on its shape. The shape parameter (the concentra- 
tion c) is very sensitive to the data at small radii, in particular 
in the cooling core region, where the mass profile is least well 
constrained, and where there could be systematic errors due 
to the PSF/projection correction. As a result we expect, as ob- 
served, increasing statistical errors on Mg as 6 decreases, and a 
corresponding increase of the raw scatter in the Mg-T relation 
(reflecting the scatter in c, see Fig. 3 of Paper I). The increase 
might be larger than that expected purely from the increase of 
statistical errors due to intrinsic scatter in c and/or systematic 
errors on c. 

4. Reliability of X-ray mass estimates 

As discussed in Paper I, there is an excellent quantitative 
agreement in shape between the X-ray mass profiles used in 
this work and the profile predicted by numerical simulations. 
The observed scaled profiles are well-described by the quasi- 
universal cusped profile (NFW-type) now found in all CDM 
simulations, and have concentration parameters as expected for 
their mass. As concluded in Paper I, this suggests that the Dark 
Matter collapse is well understood, at least down to the clus- 
ter scale. In turn, this gives us confidence in the XMM-Newton 
mass estimates, not only in the observed radial range, but also 
where we have extrapolated beyond it (i.e at 6 < <5 b s ~ 1000). 
By using the best fitting NFW model to estimate M500 and 
M2001 we nave implicitly assumed that this model remains valid 
beyond S b s (Sec. l2.lt . It would be surprising if this were not 
the case since i) it is consistent with the theoretical predictions 
above 6 b s , and ii) in one case (A1413), we were even able to 
check the validity of the NFW profile down to 5 < 500. 

Strictly speaking, the good agreement between the ob- 
served and predicted shape of the mass profiles does not mean 
that the absolute value of the X-ray mass is correct. It could be 
subject to systematic errors. However, for the correct universal 
shape of the mass profile to be recovered, this systematic error 
would have to be the same, within the statistical errors at all 
observed 6, whatever the cluster temperature. 

One possible source of such systematic error is a departure 
from hydr o static e quilibrium (HE). The recent simulations of 
iKav et all |2004a) suggest that the mass determined from the 
HE equation underestimates the true mass, due to residual gas 
motion. The effect is about the same at all radii up to 6 = 500. 
It is of the order of 15 per cent for adiabatic models and of 10 
per cent for models including cooling and feedback, with typ- 
ical variations of +5 per cent. Such variations would not sig- 
nificantly change the shape of the X-ray mass profiles, taking 
into acccount our statistical errors. Thus the measured M2500 
and Mi 000 values, and thus the corresponding normalisation of 
theM25oo-r and Miqoq—T relations, could well be ~ 10-15 per 
cent too low. The offset would be the same for the M^-T and 
M200—T relations, since they are derived from 'extrapolation' 
of the NFW model. Note that this is probably an upper limit, 
since we focus on particularly relaxed clusters. 

Another possible source of systematic error is that associ- 
ated with errors in estimated temperatures from uncertainties in 
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the instrument calibration. Possible errors are of the order of 10 
per cent. This value is consistent with the systematic difference 
observed between temperature derived with XMM-Newton and 
Chandra with the former being on a verage 0.92 ± 0.08 times 
the latter llKotov & Vikhlininl l2005). Since the mass derived 
from the HE equation scales as M oc T, error in T would trans- 
late into a systematic error on the estimate of the 'real' mass 
of clusters. However that would not change the normalisation 
of the 'observed' M-T relation, since that depends only on the 
shape of the temperature profile. 



5. Comparison with previous determinations 

In the present study, as discussed above, an NFW profile 
has been used to describe the integrated mass profile, de- 
rived, assuming HE, from the observed density and tempera- 
ture profi les. A sim ilar approach 5 was used in the Chandra 
study of lAllen et ail d200lh . Previous ROSAT/ASCA stud- 
ies also estimated the mass from the HE equation, but as- 
sumed a /?-model for the gas density profi le and a polytropic 
(or even i sothermal) temperatur e profile i Horner et al. 1999; 
Nevalainen et all200d:lFinoguenov et al. l 200lUXuet alhOOlt 
Castillo-Morales & Schindlerl 120031: ISanderson et all 12003 ) . 
Data were of poorer spatial resolution and statistical quality, 
and had less radial extent, thus requiring more extrapolation, 
particularly for low mass clusters. This could introduce sys- 
tematic errors and biases, particularly at low 6, since a mass 
profile derived from an isothermal or polytropic y6-model is 
not c onsistent with an NFW p rofile at large radii (see Fig 1 1 
of iNeumann & Arnaudlll999l) . On the oth er hand, the latest 
ROSAT/ASCA studies of the M-T relation jFinoguenov et alJ 
l200lUSanderson et alJ200l are superior in terms of the size of 
the cluster samples, and their wide and homogeneous coverage 
in temperature. Their results are compared to ours in Tabled 

Fino guenov et alJ J200ll) established the M^ao-T relation 
for 39 clusters with ASCA temperature profiles. Interestingly, 
their results are consistent with ours within the uncertainties 
(Table QJ. This suggests that systematic errors are not dom- 
inant over statistical errors at 5 — 500 in this ROSAT/ASCA 
study. However, our normalisation is on the upper side of their 
allowed values. The slope they find for their hot cluster sub- 
sample (T > 3 keV) is, as we h ave found, consist e nt wit h the 
expected a — 1.5 value. When Finoguenov et alJ J200lh in- 
cluded all clusters (down to T ~ 0.9 keV), they found a steep- 
ening of the MsM-T relation: a — 1 .78+0. 1 . The effect is larger 
than in our case: we find a — 1.71 + 0.09 (although the differ- 
ence is not significant). However, our sample does not reach 
quite such low temperatures and the difference could also re- 
flect a progressive steepening of the M-T relation toward low 
masses. 

Th e same remark holds for the results of San derson et alJ 
J2003I) . who derived a slope of a = 1.84 ± 0.06 for the M 1Q0 - 
T relation, which is barely consistent with our value. Their 



5 The method differs lightly from ours. I Allen et alJ {zOOl) predict 
the temperature profile corresponding to a given NFW mass profile 
and the observed surface brigthness profile and fit it to the observed 
temperature profile. 



Table 3. Comparison of the present results with M-T relations 
from the literature, a is the logarithmic slope of the relation and 
A is the normalisation at kT — 5 keV, in units of 10 14 M for 
H = 70 km/s/Mpc. 



Reference" 



Method b 



2500 



Observation 
Present work 
Present work 
ASF01 
SPF03 



Theory 



EMN96 
BMS04 
KSA04 



1.69 


±0.05 


1.70 ±0.07 


M NFW , T > 2.0 keV 


1.79 


±0.06 


1.51 ±0.11 


M NFW , T > 3.5 keV 


1.88 


± 0.34 c 


1.51 ±0.27 


M NFW , T > 5.5 keV 


1.2 + 


0.4 C 


1.84 ±0.14 


Mpy, T > 5.5 keV 


2.5 




1.5 


r cra , adiabatic simul. 


1.73 


±0.35 


1.55 ±0.05 


T m ,T > 2.0 keV 


1.97 


±0.03 


1.54 ±0.05 


T 

1 m 



5 = 500 



Observation 
Present work 
Present work 
FRB01 
FRB01 



Theory 



EMN96 
VBB02 
BMS04 
RMB05 
RMB05 



3.84 ±0.14 


1.71 ±0.09 


M NFW , T > 2.0 keV 


4.10±0.19 


1.49 ±0.15 


M NFW , T > 3.5 keV 


3.26 ± 0.60 c 


1.48 ±0.11 


Mpy, T > 3 keV 


3.31 ±0.45 c 


1.78 ±0.10 


Mpy, T > 0.9 keV 


5.6 


1.5 


r cw , adiabatic simul. 


3.6 


~ 1.7 


T 

1 em 


4.6 ±0.2 


1.59 ±0.05 


r era ,r > 0.7 keV 


7.2 ±0.5 


1.66 ±0.09 


T sl , T > 1 keV 


4.2 ±0.2 


1.53 ±0.05 


M Py , T A ; T > 1 keV 



200 



Observation 
Present work 
Present work 
SPF03 



Theory 



EMN96 
MTK02 



5.34 ± 0.22 
5.74 ±0.30 
4.5 ± 0.3 C 

7.4 
7.6 



1.72 ±0.10 
1.49 ±0.17 
1.84 ± 0.06 

1.5 
1.61 



M NFW , T > 2 keV 
M NFW , T > 3.5keV 
Mp y , T > 0.6 keV 

r em , adiabatic simul. 



Notes: (a): R e ference s: (ASF01) lAllen et alJ 1200 ll) : (BMS04) 
i Borgani et all <2004l): ( EMN96) lEvrard et alJ 71996); (FRB01) 
Finoguenov et alJ 1200 lh : (KSA04) iKav et alJ l2004rj); (MTK02) 
Mua nwong et alJ l2002h: fPW) present work ; (R MB05)lRasia et all 
120051) : (SPF03) ISanderson et all 120031) : rVBB02) Rfoit et alJ 120021) : 
(b): Method: (Mn F w): mass estimated using an NFW model to de- 
scribe the mass profile; (Mp y ): mass estimated using a polytropic 
/?-model for the gas distribution; (T m ): mass-weighted temperature; 
(r cm ): emission-weighted temperature; (r s i): spectroscopic-like tem- 
perature as defined in lMazzotta et alJ 120041) . (c): The normalisation 
at 5 keV is derived from the published normalisation at 1 keV and 
the best fitting slope. The fractional error has been assumed to be the 
same, which is conservative. 



large sample includes 66 clusters in the [0.5 - 15] keV tem- 
perature range. Furthermore, their normalisation at 6 = 200 is 
significantly lower than ours (by 15 per cent, Table [3}- This 
may reflect the introduction of systematic errors when extrap- 
olating polytropic models down to 6 as low as 200. However, 
such systematic errors could not explain the discrepancy be- 
tween our results and theirs at 6 — 2500. At that density con- 
trast, their slope for a hot cluster (T > 5.5 keV) sub-sample 
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is a — 1.84 + 0.14, a result which is inconsistent with ours 
at the ~ 95 per cent confidence level. Their normalisation at 
5 = 2500, (1.2 + 0.4) x 10 14 M , is lower than our value of 
(1.69 + 0.05) x 10 14 M Q and is only barely consistent with the 
Chandra results. 

ISanderson et alJ J2003I) have suggested that the discrep- 
ancy (with respect to the Chandra results) might be related 
to the dynamical state of the clusters in the different sam- 
ples. Both the Chandra study and the present XMM-Newton 
study focus on pa rticularly relaxed clusters, which is not 
the case in ISanc erson et all s study. However, we note that 
iFinoguenov et al. 's sample does not discriminate in terms of 
dynamical state, and their results are in good agreement with 
ours. That said, the M-T relation could well depend on the ex- 
act dynamical states of the clusters in the sample in question, an 
e ffect which is n ot trivial to predict. The numerical simulations 
of iRowlev et alJ (|2004) show that clusters with substructure 
tend to lie below the mean T-M relation, probably due to in- 
complete thermalisation (their Fig. 15). However, for the same 
reason, we would expect the X-ray mass to underestimate the 
true mass jKav etalj E004a). moving unrelaxed clusters back 
closer to the mean relation. A XMM-Newton study of the M-T 
relation for an unbiased sample of clusters is needed to assess 
the effect of cluster dynamical state on the measured M-T re- 
lation. The relation should ideally be compared to numerical 
simulatio ns and lensi ng mass data. In any case, we do not con- 
firm Sanderson et all s results, at least for the relaxed clusters 
considered here. 

Final ly, it is of in te rest to compare the present results with 
those of Etto ri et all ll2002l) . who also use an NFW model 
to estimate masses 6 . Their M-T relation for a relaxed sub- 
sample of 12 clusters with T > 3 keV can be directly com- 
pared with our relation for hot clusters. At 5 = 2500, they found 
a - 1.88 + 0.27, which is marginally consistent with our value, 
and a = 2.3 ± 0.4 at 6 = 500, a value clearly rejected by our 
data. However, there is a large scatter in their M-T relation at 
5 = 2500, which becomes dramatic at 5 = 500 (their Figure 
5). We do not observe such a scatter. It may reflect systematic 
errors connected to the extrapolation of the NFW model. As 
we discussed above in Sect. 13.21 the precision on extrapolated 
mass depends on the precision on the concentration parameter 
(i.e., the shape of the mass profile, especially in the center), 
which is more difficult to constrain with BeppoSAX than with 
XMM-Newton, in particular due to the larger BeppoSAX PSF 



6. Comparison with theoretical predictions 

The temperature structure of the ICM in a cluster is the result 
of the complex interplay between gravitational processes (i.e., 
the evolution of the gas in the Dark Matter potential), and of 
any other process that can affect the gas entropy (e.g., radiative 
cooling and heating from galaxy feedback). The theoretical M- 
T relation - which should be viewed rather as a T versus M re- 
lation when predicted from theoretical studies - depends on the 
exact modelling of all these processes. Moreover, as we discuss 
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6 In their approach, the predicted temperature profile was fitted to 
the observed BeppoSAX temperature profile. 



Fig. 4. The normalisation of the Ms-T relation at T — 5 keV 
for various 6. The results of theoretical works (see Table[3]for 
references) are compared to the values derived in the present 
work (PW) . 



below, comparison of observations and theory also depends on 
the exact definition of the 'average' temperature, since the gas 
is never perfectly isothermal. 

6.1. The normalisation of the M s -T relation at 5 keV 

The normalisation of the M-T relation is particularly well con- 
strained by our study, the statistical error now being less than 
~ 5 per cent at 5 keV. The value of the normalisation depends 
on the (sub)sample considered because it is correlated with the 
slope. However by choosing a reference temperature of 5 keV, 
close to the median temperature, we minimize the effect, and 
the difference is of the order of the statistical error. We can 
thus first compare our results with the predicted values, quasi- 
independently of the slope issue. 

6.1 .1 . Adiabatic models 

The difference with the prediction from adiabatic models is 
not dramatic, particularly when the dispersion among vari- 
ous adiabatic simulations is taken into account. The normal- 
isation is 30 per cent below the prediction of lEvrard et alJ 
(1996), compared to typically more than 50 per cent differ- 
ence in the normali sation derived from different adiabatic sim- 
ulations (Table 3: iHenrvl |2004). However, as discussed by 
iMuanwong et alJ d2002l) . higher resolution simulations tend to 
predict higher normalisation, thus exacerbating the discrepancy 
with our results. Adiabatic simulations probably fundamentally 
fail to predict the correct normalisation of the M-T relation. 
The observed discrepancy could in principle be due to incorrect 
modelling of the dark matter component itself, since it is this 
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which drives the potential. However, we have already argued 
that this is not the case from the excellent agreement between 
the observed and simulated mass profiles (Paper I and Sec.@}. 
A more likely explanation is that non-gravitational processes 
affect the M—T relation, as they affect all other scaling rela- 
tions. 

6.1.2. Non-adiabatic models 

From theoretical arguments and from observations of the gas 
entropy, it is now clear that both cooling and galaxy feedback 
have to be tak en into account when discussing relations inv olv- 
ingthelCM ( I Voit & Ponmanl2003l:lp"ratt & Arnaudl2005t) . We 
thus focus on published theoretical studies that include both 
phenomena. Their predictions are compared with our results in 
Tableland in Fig.0] All of the quoted simulations are made in 
the 'concordance' ACDM cosmology, and use the same defini- 
tion of Mg as in the present work. However, they do not always 
use the same definition of the temperature (see below). 

At a given mass, cooling and galaxy feedback increase the 
gas entropy as compared to the value attained by pure gravita- 
tional heating. In many scenarios, these processes increase the 
temperature and, as expected, a lower normalisation is found 
when the se processes are included in a given n umerical sim- 
ulation (Mua nwonget al.ll2002l IThomas et ai1l2002l) or ana- 
lytical model (Vo it et alJl2002l) . Recent non-adiabatic simula- 
tions seem to be quite successful at reproducin g the observed 
M*>sm -T relation. The normalisation derived bv IBorgani et alJ 
(2004), for a sample of simulated clusters with temperatures 
T > 2 keV, is in perfect agreement with our value (s ee also 
Fig. while the normalisat ion derived by iKav et al J (J2004b, 
see also iThomas et alJl2002l) is only ~ 10 per cent too high. 
These works used the mass weighted temperature, T^soOt es_ 
timated within /?2500, which is directly linked to the thermal 
energy. As discussed in Sect. 12.21 our spectroscopic tempera- 
ture, T s , should be close to the mass weighted temperature in 
that region. Thus, the good agreement between the observed 
M25W)-T S and predicted M2500— I'm relations is encouraging, and 
suggests that the thermal energy content in this central region 
is roughly correctly modelled. 

Up until recently, the standard temperature definition used 
to mimic X-ray observations in numerical simulation stud- 
ies was the emission-weighted temperature r em . Using this 
temperature definition at 6 — 500, the normalisation of 
IBorgani et alJ J2004I) is no w too high b y abou t 1 2 - 20 per cent, 
the normalisation found bv lVoit et alJ Il2002|) is margin a lly too 
low and the normalisation found by Muan wong et alJ l|2002) 
at 5 = 200 is 30 per cent too high ( Fig.0). It is likely that the 
variabality of these results is linked to differences in the various 
physical models used. 

Independent of the physics, a crucial point seem s to be the 
exact definition of the temperature. Recently, Maz zotta et alJ 
(2004) introduced the spectroscopic-like temperature (T s {) in 
order to better reproduce the temperature obtained f rom spec- 
tral fits when the ICM is multi-temperature. Maz zotta et alJ 
show that T s \ is biased towards the lower values of the domi- 
nant thermal component, and that in general T em overestimates 



r s i. Unfortunately, this exacerbates the disagreement between 
observed a n d sim ulated M500— T normalisations. Using r s i, 
iRasia et all J2005I) over-predict a normalisation , relative to our 
values, by a factor as large as ~ 1.8. Note that Rasia et alJ use 
the same physical model as Borg ani et alJ J2004I) and vet their 
normalisation is ~ 50 per cent higher (see Tableland also their 
Fig 2) . We note that the earlier work of Mathiesen & Evrard 
(2001), based on spectroscopic temperatures of adiabatic nu- 
merical simulations, showed a smaller effect. Strictly speaking 
these temperatures were estimated with R500, whereas our spec- 
troscopic temperature measurement T s is interior to 0.5i?20o(~ 
/?iooo)- A1413 is the only cluster for which we have data up to 
6 = 500. The spectroscopic temperature within 5 = 500, T St $oo 
is only slightly smaller (by 3 per cent) than T s . This would in- 
crease the normalisation of the M^qq-T relation by less than 5 
per cent if we used r S; 5oo (assuming the same correction factor 
for all clusters). 

It thus appears that there is a genuine disagreement be- 
tween observed and predicted normalisatio n of the M^ o g-T re - 
lation. One interpretation, as proposed by Rasi a et all (J2005), 
is that the X-ray mass undere stimates t he 't rue' mass (see also 
IBorgani et all2004llMuanwong et all2002l) . Using T s \ temper- 
atures, they estimated the value of M500 that an X-ray observer 
would derive from their simulation using the HE equation and 
a polytropic y6-model. The resulting M500 - T^i is indeed now 
in good agreement with our observation (Table [3] and Fig.0J. 
However, the normalisation is 4.2 ± 0.2 X 10 14 M , as com- 
pared to 7.2 + 0.5 x 10 14 M when using the 'true' theoretical 
mass. This corresponds to a very serious underestimate of the 
mass by X-ray observations: the 'true' M500 mass of clusters 
would be a factor ~ 7.2/4.2 = 1.7 higher than the X-ray mass. 
We think this is very unlikely, at least for the masses estimated 
as in the present work. Firstly, our approach - fitting an NFW 
model and extrapolating the mass profiles - is more sophisti- 
cated than the simple polytropic j6-model approach. Secondly, 
we note again the excellent quantitative agreement of our mass 
profiles with theoretical predictions. If we have underestimated 
the 'true' M500 by a factor of 1.7, we should also have under- 
estimated the 'true' M2500 by the sam e factor. Thi s is unl ikely: 
from combined lensing/X-ray studies, Alle n et alJ (12001) con- 
clude that systematic uncertainties are less than 20 per cent 7 . 
Conversely, for the X-ray mass profiles to have the correct uni- 
versal shape, as we have observed, the predicted difference be- 
tween the X-ray mass estimates and the true mass should be 
roughly constant with radius. This is not what is expected if 
the difference is important at 6 — 500: it is linked to differences 
between the temperature profile derived from projected T s i(r) 
values and the true profile, which depends on the ICM struc- 
ture a long the line of sight ( Ma zzotta et alJl2004c IRasia et alJ 
2005), and thus a priori on radius. It would be interesting to 
check this point with numerical simulations. 

A more likely explanation is that numerical simulations do 
not correctly describe the gas thermal structure at large scale, at 
least for relaxed clusters considered here. We note that numer- 
ical simulations predict temperature profiles decreasing with 

7 Note that our M2500-T relation is th e same, both in slope and nor- 
malisation, as that found bv lAllen et alJ 1200 ll) . 



10 



Arnaud et al.: The M-T relation 



radius, by nearly a factor 2 at 5 — 500 ( Bor gani et al.ll2.004l 
Fig 6), while the observed profile is flatter (Fig.0. This would 
bias low the T s \ as compared to the X-ray temperature of real 
clusters, and thus increase the normalisation of the M^qq-T^ 
relation 8 . It would be interesting to compare the theoretical 
M2500-7si and A^soo-^m relations, and investigate if there is 
continued good agreement with observations. We expect this 
to be the case since the predicted temperature variations are 
not dramati c - less than 20 per cent variations within 5 — 2500 
(see Fig 6 of Borgani et al. 2004 and Fig 9 of Kav et al. 2004a), 
so that T s i should be close to T m . 

6.2. The slope of the M-T relation 

The observed M-T relation slope is consistent with the self- 
similar expectation for the sub-sample of hot clusters (T > 
3.5 keV): or = 1.51 ± 0.11 at 6 = 2500, where it is best con- 
strained. The slope is significantly higher when the whole sam- 
ple (T > 2 keV) is considered: a = 1.71 ± 0.07. 

A value of a = 1.5 is expected from the virial theorem 
if clusters obey self- similarity. All adiabatic simulations con- 
firm this value fevrard et al Jl 1 994 IPenll 1 998t lEke et alJll998t 
Brvan & Normanlll 998l lYoshikawa et al.ll200(i iThomas et al J 
2001), including when a wide bandpass spectral tempera- 
ture, as measured with Chandra or XMM-Ne wton, is used 
to establish the M-T relation ( Mat hiesen & EvrardboOll) . 
Numerical simulations including cooling and feedback do pre- 
dict a slightly higher slope. However, the effect is smaller than 
we observe (A(a) = 0.05 - 0.1) and is generally not significant 
(Table Ell. The only ex ception is the M5oo-r s i relation derived 
bv lRasia et al.l ( 120051) : a = 1.66 + 0.09; however the normal- 
isation is then much too high (see above). It is also worth not- 
ing th at the phenomenological analytical model of IVoit et alJ 
(2002) yields a steeper slope. We obtained a ~ 1.7 by fitting 
their relation (their Fig 22) in our temperature range. This is in 
good agreement with the observed value; however, in this case 
the larger slope is mostly due to the variation of the c oncentra- 
tion of the Dark Matter with mass in their model JVoit et alJ 
l2002h . which is larger than we observe (Paper I). 

As a final remark, we want to emphasise that the observed 
discrepancy with the standard self-similar value is actually 
small. The slope increase, observed when including cool clus- 
ters, is significant at most at the ~ 85 per cent confidence level. 
Furthermore, at 2 keV, the limiting temperature of our sam- 
ple, this corresponds to only -20 per cent difference in mass 
as compared to the extrapolation of the best fitting M-T rela- 
tion for hotter clusters (see also Fig.[3j. There is scatter in the 
M-T relation, and our sample comprises only 4 cool clusters. 
We thus cannot exclude that the steepening is an artefact of our 
particular choice of clusters. We also note that the quality of 
the power law fit decreases when including low mass systems. 
This may indicate that the M-T relation is actually convex in 
the log-log plane, either across the entire temperature range, or 
below a 'break' temperature. We lack clusters in the interme- 
diate temperature range to assess this issue. Clearly, a possible 

8 Since T cm is less affected, the predicted Msoo-T sm relations would 
be in better agreement with observations. 



discrepancy between predicted and observed slopes needs to be 
confirmed and better specified by considering a larger cluster 
sample. 

7. Conclusion 

Using a sample of ten relaxed galaxy clusters observed with 
XMM-Newton, we have calibrated the local Mg-T relation, in 
the temperature range [2 - 9] keV, at four density contrasts, 
6 = 2500, 1000,500,200. We used the spectroscopic tempera- 
ture estimated within 0.5/?200 (S ~ 1000), excluding the cooling 
core region, and derived the masses at various 6 from NFW 
profile fits to precise mass profiles measured up to at least 
6 = 1000. We argue that our measured masses are particularly 
reliable. The logarithmic slope of the Mg-T relation is the same 
at all 6, reflecting the self-similarity of the mass profiles. The 
slope is well constrained and is consistent with the standard 
self-similar expectation, a - 1.5, for the sub-sample of hot 
clusters (T > 3.5 keV). The relation steepens to a ~ 1.7 when 
the whole sample (T > 2 keV) is considered. The normalisa- 
tion of the M-T relation is measured with a precision better 
than +5 per cent and is 30 per cent belo w the value p r edicte d 
by the adiabatic numerical simulations of lEvrard et al.1 {T996). 

Models that take into account radiative cooling and galaxy 
feedback are now in good agreement with the observed M2500- 
T relation. We argue that remaining discrepancies at 6 = 500 
and lower are more likely to be due deficiencies in models of 
the ICM thermal structure, to which the spectroscopic-like tem- 
perature seems to be very sensitive, rather than to an incorrect 
estimate of the mass from X-ray data. 

More detailed comparisons are needed to understand the 
origin of the discrepancies between the predicted and observed 
M-T relations. Our directly measured Miqoq-T relation now 
provides the most direct constraint at large scale for numer- 
ical simulations. Simulations of mass profiles, as would be 
determined by an X-ray observer using modern Chandra and 
XMM-Newton techniques, are also needed. This would be par- 
ticularly interesting for relaxed cluster sub-samples, and using 
better representati ons of observ ed temperatures (e.g. as pro- 
posed by iMazzotta et al.lEo04l) . Such data could be directly 
compared to observed mass profiles. This would provide infor- 
mation on i) possible overall systematic errors in X-ray mass 
estimates, and ii) further test the reliability of simulations to 
correctly reproduce the ICM structure. 

On the observational side, study of a much larger, unbiased, 
sample is needed to i) determine the exact shape of the local 
M-T relation; ii) study its intrinsic scatter, and iii) assess the 
effect of cluster dynamical state on the M-T relation. 
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